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Abstract 

The influence of tlie electrons interaction with longitudinal acoustic phonons 
on magnetoeletric and spin-related transport effects are investigated. The 
physical system under consideration has been assumed to be a two dimen- 
sional electron gas system with both Rashba and Dreesselhaus spin-orbit 
couplings. In the works which have been previously performed in this field it 
has been shown that in the non-equilibrium regime, the Rashba and Dreessel- 
haus couplings can not be responsible for spin-current where in the absence 
of other interactions such as lattice vibrations identically vanishes. In the 
current work we have employed a semi-classical method by using the Boltz- 
mann approach. It was shown that the spin-current of the system, in general, 
does not go all the way to zero at any value of the spin-orbit couplings. It 
was also shown that the spin accumulation of the system can be influenced 
by the electron-phonon coupling. 
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spin current. 



1. Introduction 



Spintronics, has been subject to many investigations and has attracted 
more and more attention, from both theoretical and experimental sides, dur- 
ing the last several years. Effective control of spin polarized transport is very 
important specially for practical applications for example in multilayers and 
many studies have been conducted to explain this phenomena l|, l2|. 
According to the results of these studies, electron spin manipulation can be 
realized by applying magnetic fields or the Rashba interaction that arises 
from the inversion asymmetry in the system and can be effectively controlled 
by applying a gate voltage js, 0, Q . The most popular method for manipula- 
tion of electron spin is applying the Rashba spin-orbit interaction (SOI). This 
kind of spin-orbit interaction plays a central role in Datta and Das spin field- 



effect transistor (SFET) 



61]. Generally SOI has a significant role in magne- 



n 



toresistance effects known as weak localization [7]. There are many interest- 
ing features that have been demonstrated for this type of spin-orbit coupling 



(SOC) in the field of spin-transport [o], |10|. For example it has been ver 



ified that spin-orbit scattering can induce localization/anti 



11 



oca 



ization tran- 



12j. Meanwhile 



sition in a two-dimensional electron gas (2DEG) system 

Rashba interaction has been also suggested for spin interference devices and 


spin-filters [13|, 



14 



The other spin-orbit coupling that provides us a new spin-dependent pa- 
rameter which should be considered to design spin- dependent devices, is the 
Dreesselhaus coupling. As it is known, Dreesselhaus coupling is induced by 



2 



the bulk inversion asymmetry. The effects of spin-orbit couphngs (SOCs) in 
semiconductors have attracted growing interest due to their rules in semi- 
conductor spintronics. 

Manipulation of spin makes new functionality in electronic devices. Control 



of spin accumulation by 



field of spintronics 



16 



spm-or 



17 



18 



3it interactions has a great potential in the 



19|. In the presence of these two different 



spin-orbit interactions, i.e. the Rashba and Dreesselhaus couplings, in a two 
dimensional electron gas system, one can effectively control both magnitude 



18j. Meanwhile spin cur- 



and direction of non-equilibrium spin accumulation 
rent vanishes exactly, which implies that no spin-polarized current accom- 
pany with the spin accumulation of the system in non-equilibrium regime 
induced by an in-plane driving electric field 18 1. 



Rashba obtained non-vanishing spin-current in equilibrium state, therefore 
this spin-current can not describe any real transport of spins in non-equilibrium 
regime induced by an in-plane driving electric field |20j . Meanwhile for a non- 
equilibrium system highly anisotropic spin response to an in-plane electric 



field has been discovered 



2l|]. However, as mentioned before based on the 



semiclassical approach, Huang and Hu have shown that the non-equilibrium 
spin current vanishes exactly in two-dimensional electron gases, in the pres- 
ence of both Rashba and Dresselhaus couplings, although one can obtain 
non-equilibrium spin accumulation in this case 18(|. In addition Inoue and et 
al obtained similar results based on the Green's function approach for two- 
dimensional electron gases 22 1 . 

n 

It should be noted that the work describe by Huang and Hu [18|, has been 
based on a semi-classical Boltzmann approach developed by Schliemann and 
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Loss in [23|. They formulated the anisotropic effects of the energy disper- 
sion relation and scattering matrixes in the presence of spin-orbit couplings. 
Meanwhile the exact solution to the Boltzmann equation for two-dimensional 



181 ] in the same 



anisotropic systems has been performed by Vyborny and et al [2J], where it 
was shown that for a Rashba type two-band model the discrepancy between 
the exact and approximative Schliemann and Loss approach, remains only 
on the level of the higher order corrections [24]. Therefore we have employed 
the Schliemann and Loss approach so that the results of the current work 
can be comparable with the results of the Huang and Hu 
theoretical framework. It can be easily shown that in the presence of both 
of the spin-orbit couplings i.e. Rashba and Dresselhaus interactions, Schlie- 
mann and Loss method is still a good approximation in comparison with the 
exact Vyborny approach. This is due to the fact that the anisotropic term 
of the energy dispersion relation is negligible and the anisotropic effects can 
be inter onely through the scattering matrixes. 

In this work, based on the mentioned semiclassical approach, we have consid- 
ered the influence of electron-phonon scattering on spin-transport character- 
istics of a two-dimensional electron gas with forgoing spin-orbit interactions. 
We have verified that the electron-phonon scattering results in non-vanishing 
spin-current that can be controlled by spin-orbit interactions. It has also 
been found that the spin-current is influenced by the electron-phonon cou- 
pling. The details of the numerical results have been briefly addressed in the 
present work. 
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2. Model and approach 

The total Hamiltonian is given by 

H = Ho + Vim + Hel-ph, (1) 

in the above Hamihonian, Hq is the kinetic energy and spin-orbit interactions 
(including both Rashba and Dresselhaus spin-orbit couplings), for a 2DEG 
namely 

Ho = + a(aa:ky - ayk^) + ^{cr^k^c - ^yky)^ (2) 

where k is the wave vector of conduction electrons, ai{i = x,y) are Pauli 
matrices and a and f3 denote the Rashba and Dresselhaus strengths, respec- 
tively. 

For a given wave vector k, 

are the eigenfunctions of Hq where A = ±1 and 0k is defined easily by 

tan0 — '^^^ 

aky + /3kx 

The corresponding eigenvalues of are 

ekA = ^ + X^{a^ + l3^)k^ + AaPk,ky. (5) 

The expectation values of the electrons spins along to the x and y directions 
in the state | kA > cab be easily find to be 

'5S(k) = ^Acos(0k), 5i^(k) = ^Asin(0k). (6) 

5 



The second term of Hamiltonian, Vim(r) is the impurity scattering potential, 
and can be defined by 

VUr) = Y.^J^-rfi{vmv-v^, (7) 

i 

where the sum is performed over all of the randomly distributed impurities, J 
is the exchange interaction strength of magnetic impurities with conduction 
electrons and m(r) is the unit vector along the local magnetization. 

/ I J 3 J X ^ J y \ 

< k A I Vira{r) \ kAn^ >= Ck'.kf^n^n, . (8) 

In which 

Jz = m^J, Jy = niyj, = m^J, (9) 

and 

Ck',k = {l/^/LJ7y)J2^M^ik' - k).fj). (10) 

j 

For long range magnetic interactions, because of the shape anisotropy, we 
take rriz = and for randomly oriented magnetic moments of impurities, one 
can assume 

< rrix >=< rriy >= 0, < ml >=< rriy >= -. (11) 

Since it was assumed that the magnetic moments of the impurities have 
been randomly oriented therefore they can not be responsible for any spin- 
polarized effect such as spin current. 

The last term of Hamiltonian, Hei^ph is the electron-phonon interaction and 



can be expressed as 



25|, 



Hel-ph = ^acV.u(r), (12) 



here, Dac is defined as deformation potential for electron scattering by acous- 
tic phonons and u{r) is the small displacement vector of an ion from its 
equilibrium position, R. 

for the two-dimensional system, the displacement is determined as 



u(r 



where, M and N are mass and number of ions, respectively. Sq is the unit 
vector in displacement direction and Uq = VsQ in which Vg is the sound 
velocity and q is the wave vector of the longitudinal acoustic phonons. 
Using from eq. (??), the electron- phonon interaction can be written as 



by Defining 



we obtain the following result, directly 

Hei-ph = J][c(g)a,e^^"+ c*(g)ate-'^^1. (16) 
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The eigenstate of the phonon Hamiltonian in harmonic approximation is 
defined by \nq >, where Uq = ^ , is the Bose- Einstein distribution 
function of phonon, so we can define a new basis as follows, | kAn^ >=| 
]s.X > <S> \ riq > . Scattering matrix of electron-phonon interaction is given as 
follows 

Sn'^^ng-i6x'\c{q)./n^, if k' = k + q, 

< k'X'n'q I Hei-ph I kAng > = I (17) 

5n'^,n^+lSx'xc*{q)^/7^ + T, if k' = k - q. 

7 



n this work we have employed the procedure which have been used in 



18|. The two last terms of the Hamiltonian are responsible for both spin- 



dependent and spin-independent relaxation mechanisms. If we rename these 
two terms as 

V = Vim + Hel-ph (18) 

The Lippman- Schwinger scattering state of a conduction electron reads, 



, , I , , ^'A'n' kAnq . , , 

kXug >scat = I kArig > + 2^ _ | k'g'A' > 



(19) 



where rj is a small positive quantity. Then the expectation value of the 
electron spin in a given scattering state is 

kw ekA-ek'A'+^^ 

Here we have defined 5'A,j(k) =< kXrig \ Si \ kXrig >scat- Therefore one can 
obtain spin expectation value as seen in eq. (12T]). 



Sx,i{k) = S'-^lik) + ^ Re[Vk'\'ng,k\nq < ai >k'A'n„kAn, Pr ^ 

— Vk'X'ftgM^ng < >k'X'nq,kXnq ^T^^i^kX — Cfc'A')]- (21) 



Here < cTj >k'A',kA is the expectation value of the Pauli Matrix for a Lippman- 
Schwinger scattering state. Then, net spin density can be giving by 

<5,>=5^5A,(k)/A(k,q), (22) 

k,gA 

in which /^(k, q) is the non-equilibrium distribution function of conduction 
electrons. In the absence of external electric field, this is reduced to the 
equilibrium Fermi-Dirac distribution, 

/,(k, q) = /o(ekA) = 1^^. (23) 

l + e ^'bt 
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We have used the Debye model, so the summation over q is easily calcu- 
lated by replacing it with an integral. This integral can be considered to be 
evaluated in the interval starting from q = to the Debye wave vector, q/5. 
This wave vector is directly related to the free-electron Fermi wave vector. 
In two-dimensional metals, qD = -y/f ^f? where kp is the free-electron Fermi 
wave vector, and z is the nominal valence |26|. 

In the presence of scatterings, the non-equilibrium distribution function will 
be derived by solving the Boltzmann equation (in steady state for a homo- 
geneous system), 

^ ,f )„„. (24, 

where k = and {^)coii is called the collision integral due to scatterings, 
that in elastic scattering approximation reads |27l] 
d fx 

k'q'A' 

A'nq,kAnq /v(k',g')(l - /A(k,g))5(ekA - eyy).{25) 

k'g'A' 

Wk'A'ng.kAn, are the transition probabilities that are given by the Fermi's 

golden rule, Wk'X'nqMn, = xl^k'A'ng.kAnJ^- 

since ^n^,n, selects only the diagonal elements of |Vk/A'ng,kAngP while ^n^.n^-if^k'.k+q 
and Sn'^^nq+i^k' select some of the non-diagonal elements of | Vk'A'ng.kAn^P, 
therefore one can easily obtain 

l^k'A'nq.kAnJ^ = | < k'AV, | Hel-ph \ kAn^ > ^ + | < k'X'u'g I Virr, \ kXUg > ^,(26) 

and accordingly 

W^k'A'ng,kAnq = ^k'A'ng.kAn.g + ^k'A'ng,kAnq + ^k'A'ng,kAn, • (^7) 

9 



^k'A'nq.kAn, IS due to impurity potential, 



<l'n,-,kAn, = 11 \ ] n^Sn'^,n,. (28) 



In which we have used the following approximation, (l/LxLy) exp{i{k'- 
~ "0')) ~ "^ii "where Ui is the impurity density, Lx-, Ly are system di- 
mensions and it should be noted that the summations have been performed 
over the random positions of impurities. 

(2) (3) 

^k'A'ng kAn, ^^"^ ^k'A'ng kAn, electrou-phonou interaction, where 

(2) 

^k'A'ng kAn, for phouou absorbtiou contribution, 

W^S'n.,kAn„ = „,_i4'.k+q<5A',A ( "'"^ ^ | , (29) 



Cq^q 



(3) 

and W^/y^^ ^^^^ should be considered for the case of emission, 



<k^,kAn, = 5n'„n,+l5k',kVA',A I ' ' ' ' " ^ |. (30) 



c*v/^^;ti 
c^y^^^n 

Since the energy dispersion of conduction electrons becomes anisotropic in 
the presence of Rashba and Dresselhaus spin-orbit interactions, this manifest 
itself in the scattering process and one can choose an anisotropic solution to 



the Boltzmann equation as follows 181] 



5/,(k, q) = e^|^|^[akqA(E. ^r^x) + bkqA(E x e,).VkA]. (31) 

VkA = -^^kCkA, (32) 

where VkA is the velocity of conduction electrons, 6f\ = /a — /o, is a 
unit vector along the direction of normal of the two-dimensional plane, OkqA 

10 



and bkqx are two unknown coefficients that can determine self-consistently by 
using eq. fl2^ and eq. fl3T|) in which ^k^A arise due to the anisotropic nature 
of *e syste,„. Theu cue ca,> tha. pe uuknowu coefficients a.,, and 
are satisfying the following equations [18| 



'kg A 'kg A 



1, (33) 



OkgA _ &kgA ^ (34) 
(2) (1) ^' V"^^ 



''"kgA ''kgA 



From Eqs. (1551) and one easily gets 



^(1) ^(2) 

C^kqA = (T) ! ^kgA = (5) ; (35) 

1 + [%f ]2 1 + [ 

''kqA kg A 



in which r^^^ and r^^^,^ are two relaxation times, and are defined by 

EM4'A'ng,kAn„{l " COs[6'(VkA A Vk'A')]}) (36) 



'kgA k'g',A' ' 



EW'k'A'ng.kAngr^^ sin[6'(VkA A Vk'A')]) (37) 



t(2) ^ ^ IVk'A 

^kgA k'g',A' ' 



where 6'(vkA A Vk'A') is the angle between VkA and Vk'A'- The spin current 



operator is defined as 



20| 



(38) 



where Vx = h ^(|^) is the velocity operator. The expectation value of the 
spin current in a given scattering state i.e. J* (k, A) =< kXug \ J* | kXug >scat 
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can be obtained as shown in eq. (15^ . 

ji{k,X) = jTHKX) + hy^Re[V^>yn,Mn,<Jl>^'X'n,Mn,Pr 

— Vk>X'n kXn„ < Jx >lL'\'n„,k\nq ^T^^i^kX — ^k'X')]- (39) 



In which we have defined Jx^-' =< kAn^ | J* | kArig >. 

Then the transport spin current in x direction with spin parallel to the x and 
y axes, is given by 

J: = J2 -^-(k, A)5A(k, q), {t = x,y or Sx, Sy) (40) 

3. Results 

In the current work the spin accumulation and spin-current of a two di- 
mensional electron gas system have been obtained in the presence of the 
Rashba, Dreesselhaus and electron-phonon interactions. This has been ac- 
complished by utilizing a semi-classical model developed for anisotropic sys- 
tems. This anisotropy can be induced by spin-orbit couplings in scattering 
matrix or in the energy dispersion. As mentioned before the Rashba and 
Dreesselhaus couplings can not be responsible for spin-current generation in 



non-equilibrium regime [18|, |22|. Meanwhile non-equilibrium spin accumu- 
lation can be effectively controlled by these spin-orbit interactions [1^. In 
the present work it was shown that the electron-phonon interaction has a 
considerable role in the generation of spin-current which was expected to be 
obtained by spin-orbit couplings. The spin accumulation of the system has 
also been controlled by the strength of the electro-phonon coupling. 
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In the present system electric field is assumed to be applied along the x di- 
rection and the numerical parameters have been chosen as follows = lOeV 
is the Fermi energy, J = O.leV, rti = lO^^cm"^ is the density of impurities, 
T = IK and Vg = 4950m/s. In addition Rashba and Dreesselhaus couplings 
have been denoted by = ma'^/h'^ .ep = m[3'^ /h^ respectively. The Rashba 
coupling can reach high values up to 0.2 eV for example in epitaxial graphene 
grown on a Ni(lll) substrate 28|]. However in the present work a typically 
lower range has been chosen for SOC as reported for other materials. 
It was found out that the electron-phonon interaction can not be considered 
as a underestimate effect on spin-dependent mechanisms. It was demon- 
strated that at low electron-phonon coupling strengths the lattice vibrations 
are more effective. 

In Fig. [5] and Fig. [5] longitudinal spin-current has been depicted as a func- 
tion of the deformation potential. Spin-current has been induced due to the 
lattice-electron interactions. These figures clearly show that spin-current of 
the system has a accountable value in which its sign and magnitude can be 
controlled by the SOCs. At the same time, as shown in Figs. [5] and |5] longi- 
tudinal and transverse spin accumulations can be effectively changed by the 
Rashba and Dreesselhaus couplings. 

Therefore forgoing results show that when the effect of electron-phonon in- 
teraction is taken into account, in the semi-classical regime it turns out that 
the both components of the spin-current can take non-zero values. It seems 
that the details of the scattering potential has an important role in the the 
generation of the spin current in the presence of the Rashba and Dreessel- 
haus couplings. As reported in |l8| Huang and Hu found that short-range 



13 



delta function impurity scatterings (where actually have spherical symme- 
try) result in zero spin-polarized current in the system. As mentioned before 
anisotropic effects can be induced by two different sources: the energy dis- 
persion relation and the scattering matrix of the relaxations. The first source 
of the anisotropic is small and can be neglected. Meanwhile the anisotropic 
effects of the scattering matrix can be regarded as anisotropy of the eigen- 
states and the anisotropy of the scattering potential. For isotropic scat- 
tering potentials the anisotropy of the response functions are just given by 
the eigen-states. Therefore redistribution of the carriers' population by the 
anisotropic or inelastic relaxation mechanisms can change the ensemble aver- 
age of Rashba and Dreesselhaus k-dependent effective field. This can produce 
non- vanishing spin current in the system since the effective field in the present 
case has been modified by the in elastic scatterings. Increasing the electron- 
phonon coupling strength decreases both spin current and spin accumulation 
of the system as shown in Fig. O- Fig. El Unlike the intermediate range of 
the deformation potential at high electron-phonon couplings the momentum 
of the electrons has been effectively randomized by the electron-phonon in- 
teraction since the relaxation time of the states decreases. Therefore in this 
case, the population of the carriers approaches to the limit of the isotropic 
scatterings in which the spin current of the system vanishes. 
An important feature which can be inferred from the results is the fact that, 
the absolute value of spin-current and spin accumulation decreases for high 
electron-phonon couplings as depicted in Fig. [5] and Fig. |5l Spin-current in- 
duced by the lattice longitudinal vibrations disappears in the limit of high de- 
formation potential and rapidly increases for low electron-phonon couplings. 
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Unlike the spin-orbit couplings, electron-plionon interaction can change the 
order of magnitude of the spin-current. However it should be noted that 
numerical results reveals that, in the limit of Dac — ?■ spin current vanishes 
abruptly, and it was a numerical discontinuity (has not included in the fig- 
ures) at Dac = 0. Numerical results show that for Dac = and in the case 
of nonmagnetic impurities, J = spin-current identically vanishes. This is 
in agreement with the results that have been pointed out in Ref. {3| for 
identical conditions. 



4. Conclusion 

In this work a semi-classical approach have been implemented for the 
study of the magnetoelectric effects of a 2DEG system. The primary focus of 
this work is to show that the electron-phonon coupling has an important role 
for generation of the spin-current in non-equilibrium regime. Since it was ver- 
ified that the Rashba and Dreesselhaus couplings (when the electron-phonon 
coupling is absent) can not be responsible for spin-current in this regime. 
It was numerically verified that even at low electron-phonon couplings the 
lattice vibrations can mediate in the spin-transport process modulated by 
spin-orbit interactions. 
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Fig. 1: Longitudinal spin current as a function of the deformation poten- 
tial for different SO couplings. 

Fig. 2: Transverse spin current as a function of the deformation poten- 
tial for different SO couphngs. 

Fig. 3: Longitudinal spin accumulation as a function of the deformation 
potential for different SO couplings. 

Fig. 4: Transverse spin accumulation as a function of the deformation po- 
tential for different SO couplings. 
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